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Abstract 

We consider the mass spectrum of the 3d 3-state Potts model in the broken phase 
(a) near the second order Ising critical point in the temperature - magnetic field 
plane and (b) near the weakly first order transition point at zero magnetic field. In 
the case (a), we compare the mass spectrum with the prediction from universality 
of mass ratios in the 3d Ising class; in the case (b), we determine a mass ratio to 
be compared with the corresponding one in the spectrum of screening masses of the 
(3+l)d SU(3) pure gauge theory at finite temperature in the deconfined phase near 
the transition. The agreement in the comparison in the case (a) would represent 
a non-trivial test of validity of the conjecture of spectrum universality. A positive 
answer to the comparison in the case (b) would suggest the possibility to extend 
this conjecture to weakly first order phase transitions. 
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1 Introduction 



Universality is a powerful concept since it establishes the common long-distance behavior 
of theories characterized by different microscopic interactions, but possessing the same un- 
derlying global symmetry. A remarkable example is provided by (d+l)-dimensional SU(N) 
pure gauge theories at finite temperature which undergo a confinement/deconfinement 
phase transition associated with the breaking of the center of the gauge group, Z(N), the 
order parameter being the Polyakov loop [TJ [2] - When the transition is second order, the 
long-range critical behavior of (d + l)-dimensional SU(N) pure gauge theories at finite 
temperature is conjectured to be the same as the ci-dimensional Z(N) spin model in the 
critical region near the order /disorder phase transition [3]. As a consequence, the gauge 
theory and the spin model are predicted to have the same critical indices, amplitude ra- 
tios and correlation functions at criticality. This prediction has been accurately verified 
in several cases - see, for instance, Refs. [U El El El El [10] and, for a review, Ref . [11] . 

A few years ago a study of the broken symmetry phase of the 3d Z(2) (Ising) class 
has brought compelling evidence that universality has a much wider reach than usually 
expected. In particular, it has been shown that the Ising model and the lattice regularized 
(f) A theory both exhibit a rich spectrum of massive excitations and that mass ratios coincide 
in the scaling region [121 IS] • This result is quite far from obvious: since only the lowest 
mass contributes to the free energy, there is no simple reason why higher masses in the 
spectrum should be universal. Later on, numerical evidence has been given that the 
same spectrum characterizes (3+1) d SU(2) gauge theory at finite temperature in the 
scaling region above the deconfinement temperature [TH [15] . In the deconfined phase of 
the SU(2) gauge theory the counterpart of the massive excitations in the broken phase 
of the 3d Ising model are the so-called screening masses, i.e. the inverse exponential 
decay lengths of correlation functions between suitably defined operators, built from the 
Polyakov loops [HI [T5] . 

It would be quite interesting to extend the test of the universality of the spectrum 
in two different directions: (1) by considering a theory with a global symmetry different 
from Z(2), which, however, presents in its phase diagram a second order critical point in 
the 3d Ising class; (2) by verifying if and to what extent the universality of the spectrum 
holds also in the broken symmetry phase near a weak first order phase transition. 

The 3d 3-state Potts model with external magnetic field model provides a good test- 
field for both these investigations. Indeed, its phase diagram in the temperature - magnetic 
field (T, h) plane (see Fig. [T|) exhibits a line of first order phase transitions which starts 
at (0,0), moves along the T-axis up to a transition temperature T t and then bends in 
the positive /i-plane up to reaching a second order endpoint, which belongs to the Ising 
class [16J. The 3d 3-state Potts model near this critical endpoint should exhibit a mass 
spectrum in the Ising universality class. Instead, the mass spectrum of the 3d 3-state 
Potts model at zero magnetic field, in the broken phase near T t , i.e. for T < T t should 
reproduce the spectrum of screening masses of the finite temperature (3 + i)d SU(3) pure 
gauge theory in the deconfined phase near criticality, if universality holds also in the case 
of weak first order phase transitions. 
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Figure 1: Qualitative phase diagram of the 3c? 3-state Potts model: the solid line in bold 
is the line of first order phase transitions. T t is the order /disorder transition temperature 
at zero magnetic field, (T c , h c ) is the endpoint in the 3c? Ising class. 

In this paper we determine the low-lying masses of the spectrum of the 3c! 3-state Potts 
model in two different sectors of parity and orbital angular momentum, + and 2 + , and in 
two different regions of the (T, h) phase diagram, (a) near the critical Ising endpoint, (b) 
near the transition temperature T t at zero magnetic field. For the determinations in the 
region (a), the comparison between the resulting mass ratios and the corresponding ones in 
the 3c! Ising model [12] will provide a check of validity of the Ansatz of universality of the 
spectrum. The mass ratios determined in the region (b) are instead the quantities to be 
compared with the corresponding screening masses in the finite temperature SU(3) gauge 
theory in order to verify the validity of the Ansatz near a weakly first order transition 
point. 

This paper is organized as follows: in Section [2] we briefly review the 3c! 3-state Potts 
model and describe its phase diagram on the temperature - magnetic field plane; in 
Section [3] we describe the method used to extract the massive excitations in the two 
regions of the phase diagram in which we are interested; in Section H] we present our 
numerical results and in the last Section we draw our conclusions. 



The 3-c! 3-state Potts model [T71 [18] is a spin theory in which the fundamental degree 
of freedom, Sj, defined in the site i of a 3-dimensional lattice, is an element of the Z(3) 



2 The 3d 3-state Potts model 



group, i.e. 



Si = e's 



Oi = {0,1,2} 



(1) 



The Hamiltonian of the model is 
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up to an irrelevant constant. Here, (3 is the coupling in units of the temperature and the 
sum is done over all the nearest-neighbor pairs of a cubic lattice with linear size L. 
The Hamiltonian (T5]) is invariant under the Z(3) transformation 

Si — > s\ = e^Si , (3) 

where a is fixed at any of the values {0, 1, 2}. It is well known that this system undergoes 
a weakly first order phase transition [19], associated with the spontaneous breaking of the 
Z(3) symmetry. The order parameter of this transition is the magnetization, 

i 

which, in the thermodynamic limit, is zero above the transition temperature T t (or below 
the transition coupling (3 t ) and takes a non-zero value below T t (or above f3 t ). 

In presence of an external magnetic field it is convenient to work with an Hamiltonian 
written in terms of the <j\ degrees of freedom, instead of the Sj ones. For a uniform 
magnetic field along the direction with strength h in units of the temperature, the 
Hamiltonian is 

(ij) i 

where E is the internal energy and M is the magnetization. The magnetic field breaks 
explicitly the Z(3) symmetry. However, first order transitions still occur for values of the 
magnetic field strength h below a critical value h c , the transition temperature increasing 
with increasing h. The line of first order phase transitions ends in a second order critical 
point P c = (T c , h c ) (see Fig. [1]), belonging to the 3d Ising class |16j . The Hamiltonian in 
the scaling region of the critical point P c can be written as 

H = -tE - £M , (6) 

where E and M are the Ising-like energy and magnetization operators and r and £ the 
corresponding temperature-like and symmetry-breaking-like parameters. This means that 
(M) is the new order parameter. Close enough to P c , the following relations hold, 

M = M + sE , E = E + rM , (7) 

where the mixing parameters (r, s) have been determined numerically for several lat- 
tice sizes L in Ref. [16J. The r-direction identifies the first order line (see Fig. [2]). In 
Ref. [16] the location of the critical point P c has been accurately determined: P c = 
(/3 c ,/ lc ) = (0.54938(2),0.000775(10))or, equivalently, P c = (r c , £ c ) = (0.37182(2),0.25733(2)). 

3 Massive excitations and Universality 

Among the quantities relevant in the description of a phase transition there is the cor- 
relation function of the local order parameter: in the case of the 3d 3-state Potts model 
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Figure 2: Qualitative phase diagram of the 3d 3-state Potts model in the (/3, /i)-plane: £ 
and r are the symmetry-breaking and the temperature parameters of the Ising theory ([6]). 



this is just the local spin Sj. The point-point correlation function is defined as 

r i0 (r) = (s\s io )-(sl)(s io ), (8) 

where i and io are the indices of the sites and r is the distance between them. The large-r 
behavior of the point-point correlation function is determined by the correlation length of 
the theory, £ , or, equivalently, by its inverse, the fundamental mass. In order to extract 
the fundamental mass it is convenient to define the wall-wall correlation function, since 
numerical data in this case can be directly compared with exponentials in r, without any 
power prefactor. The connected wall-wall correlator in the x-direction is defined as 



G(x) = — '^2(w(xq + x) t w(x )) - (w(xq + (w(x )) 



J-t) 



where 



w(x) 



J2 S i X 'Vr 



} 



(9) 



(10) 



represents the spin average over the "wall" at the coordinate x. 
The general behavior for the function G(x) is 



G(x) = y^a w e" 



■m n x 



where mo is the fundamental mass, while mi, ... are higher masses with the same 
angular momentum and parity (0 + ) quantum numbers of the fundamental mass. Mass 
excitations in channels different from + can be determined by a suitable redefinition of 
the wall average (flQj) . The fundamental mass in a definite channel can be extracted from 
wall-wall correlators by looking for a plateau of the effective mass, 



m eS (x) = In 



G{x) 
G(x - 1) 



(12) 
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at large distances. Excited masses in the given channel can be found by the variational 
method [2U1 EE] , which consists in defining several wall-averaged operators with the given 
quantum numbers and building the matrix of cross-correlations between them. The eigen- 
values of this matrix are single exponentials of masses in the given channel. The possibility 
to determine masses above the fundamental one in a given channel relies on the ability 
to define operators which have a large overlap with the excited states. Usually, these 
operators are defined in order to probe different length scales. In the present work we 
considered only the + and the 2 + channels; the local variables to be wall-averaged as 
in (TTOl) have been defined in the following way: 

S {x,y,z}( n ) ~ s {x,y,z}( s {x,y+n,z} + s {x,y,z+n}) > 

S \x,y,z}( n ) — s {x,y,z}(S{ Xj y +riiZ } — S{ Xj y :Z+n }) . (13) 

Typically, we have used about ten operators in each channel, with the largest n almost 
reaching L. We can anticipate that we have not been able to extract masses higher than 
the fundamental one in any of the two channels. However, the use of the variational 
method has improved considerably the evaluation of the fundamental 2 + mass. Need- 
less to say that the determination of the fundamental mass in the + channel by the 
definition (fT3"|) turned out to be in perfect agreement with that from 

A few years ago it has been proposed [12] that ratios between massive excitations 
in the broken phase are universal quantities, i.e. mass ratios must be equal in theories 
where the same mechanism of symmetry breaking at the transition is active. One of the 
motivations of the present work is to check whether this statement holds also for the 
broken phase of the Ising theory (jH]) near the critical endpoint P c of the 3d 3-state Potts 
model. This would provide us with a non-trivial check of validity of the conjecture of 
universal mass spectrum. 

For the 3d Ising universality class, a definite prediction exists for the ratio between 
the fundamental + and 2 + masses, m 2 +/m + [12J. We can estimate the same ratio near 
the Ising endpoint of the 3d 3-state Potts model and verify if it is consistent with this 
prediction. The procedure to determine the fundamental masses in the two channels of 
interest is the same outlined above, with the only difference that we need to use the correct 
local order parameter to build correlators. We have defined this local variable, rhi, in such 
a way that it reconstructs the global magnetization operator M after summation over the 
whole lattice: 

m i = 8 aiah + -J26 t7iai+fi . (14) 
A 

Indeed, it is easy to see that Mj = ^i- 

The second aim of the present work is to study mass ratios also in the broken phase 
of the 3d 3-state Potts model at zero magnetic field near the transition point. It is 
known that this transition is weakly first order, a signature of this being the fact that the 
correlation length £, though finite at the transition point, becomes much larger than the 
lattice spacing. This fact could open the door to universality effects, such as the above- 
mentioned conjecture on mass ratios. If this claim is true, an important consequence 
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would be that mass ratios in the broken phase near criticality of the 3c? 3-state Potts 
model at zero magnetic field should reproduce the corresponding ratios in the phase of 
broken Z(3) symmetry of the (3+l)c? SU(3) pure gauge theory at finite temperature. This 
is the deconfined phase of the SU(3) gauge theory and the counterparts of the massive 
excitations of the spin model are the so-called screening masses, i.e. the inverse decay 
lengths of the Yukawa-like interaction potential between static color sources. 

As a probe for the check of consistency of this scenario, we can consider again the 
mass ratio m 2 +/m +. If it will turn out that its value near the transition coupling (3 t 
in the broken phase of the 3d 3-state Potts model is compatible with the corresponding 
ratio of screening masses in the deconfined phase near transition of (3+1)6? SU(3) at 
finite temperature, this could be taken as an indication that the conjecture of universal 
spectrum can be extended also to the case of weakly first order transition. 

Summarizing, we have to calculate the m 2 + / m + ratio in two regions of interest in 
the temperature - magnetic field phase diagram: (a) in the broken phase near the critical 
Ising endpoint, (b) in the broken phase near the transition point of the model in absence 
of magnetic field. In the case (a), we have to compare the result with the predicted value 
from universality in the 3c? Ising class, given in Ref. [12]; in the case (b), our determination 
is to be compared with the corresponding one in the SU(3) pure gauge theory [22J. 

4 Numerical results 

We have performed numerical Monte Carlo simulations of the 3o? 3-state Potts model using 
a cluster algorithm [23l 124"] to reduce the autocorrelation effects. In order to minimize the 
finite volume effects, we have imposed periodic boundary conditions. Data analysis has 
been done by the jackknife method applied to bins of different lengths. 

For the simulations near the critical endpoint P c (region (a)) we used lattices with 
size L = 70; near the transition coupling f3 t at zero magnetic field (region (b)) we chose 
instead L = 48. In both cases we have seen tunneling between degenerate minima near 
the transition point. This finite volume effect can spoil mass measurements in the scaling 
region and must be treated carefully. Depending on the order of the transition, tunneling 
effects show up differently and must be removed accordingly. 

4.1 Region (a): critical endpoint 

We have performed simulations on 70 3 lattices for which the mixing parameters appearing 
in © turn to be s(L = 70) = -0.689(8) and r(L = 70) = 0.690(3) [16]. 

First of all we have considered the distribution of the order parameter M in the broken 
phase of the Ising theory ([6]) near P c . Fig. [3] shows the structure, typical for a second- 
order phase transition, with two peaks corresponding to the two degenerate minima in 
the broken phase which separate while moving away from P c along the first order Ising 
critical line. This double-peak structure is the signal of tunneling. The last plot in 
Fig. [31 obtained for Ising couplings £ = £ c and r=0. 37248, shows two almost completely 



7 



separated peaks. We have decided to choose this as our working point and performed 
here simulations with statistics 200k. 

We have removed here tunneling effects by the brute force method of analyzing sepa- 
rately data belonging to each peak. In Figs. H] and [5] we show the behavior of the effective 
masses in the + and in the 2 + channels at £ = £ c and r=0. 37248, as functions of the sep- 
aration between walls, for the configurations in the "right-peak" . Similar plots have been 
obtained for the configurations in the "left-peak" . In each case, the plateau mass is taken 
as the effective mass (with its error) belonging to the plateau and having the minimal 
uncertainty. We define plateau the largest set of consecutive data points, consistent with 
each other within la. This procedure is more conservative than identifying the plateau 
mass and its error as the results of a fit with a constant on the effective masses m e ff(x), 
for large enough x. We have found 

"right-peak" (statistics 115k) 

am + = 0.0725(63) , am 2 + = 0.1981(87) , — = 2.73(36) ; (15) 

m + 

"left-peak" (statistics 85k) 

am + = 0.0714(40) , am 2 + = 0.1959(80) , — = 2.74(27) . (16) 

The uncertainty on the mass ratios has been determined by usual propagation of errors. 
The two determinations of the masses and, therefore, of the mass ratios are consistent, 
as expected. Moreover, they are compatible with the value of the 3d Ising class |12j . 
m 2 +/m + = 2.59(4). 

4.2 Region (b): transition point at zero magnetic field 

We have performed simulations on 48 3 lattices for several values of the coupling (3 in the 
broken phase of the 3d 3-state Potts model at zero magnetic field. A summary of the 
Monte Carlo simulations is presented in Table [lj 

Close enough to (3 t {L = 48)=0. 550538, determined in Ref. [19], the scatter plot of the 
complex parameter S shows the coexistence of the symmetric phase (points around (0,0) 
in the Im S - Re S plane in Fig. [6]) and of the broken phase (points around the three roots 
of the identity in Fig. [6]). Notice that the peaks in the distribution of points in Fig. [6] 
are well separated, as it must be for first order dynamics. In the thermodynamic limit 
this would occur only at the transition point; at finite volume, tunneling between broken 
minima and the symmetric phase occurs in a small region around the transition. The 
amplitude of this region decreases with the volume. For L=48 it is of the order of 10~ 4 
in [EJ. 

Moving away from (3 t the symmetric phase becomes less and less important, up to dis- 
appearing. For /3=0.5508 there are only the three broken minima, as shown in Fig.[7](top). 
For larger j3 values the tunneling between broken minima survives. However, the three 
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peaks are totally separated and it is therefore possible to "rotate" unambiguously all of 
them to the real sector (see Fig. [T^bottom)). Working only in one sector allows us to 
optimize statistics. With this approach we have calculated the fundamental masses in 
the + and 2 + channels for several j3 values, up to 0.60 ((T t — T)/T t ~ 0.08). In Figs. [8] 
and [9] we show the behavior of the effective masses in the + and in the 2 + channel at 
(3 = 0.554, as functions of the separation between walls. Similar plots have been obtained 
for the other (3 values. The plateau mass values have been determined as described in 
the previous Subsection. In Table [1] we present our results, whereas in Fig. [10] we have 
plotted the behavior of m + and m 2 + versus f3. 

Table 1: Fundamental masses in lattice units in the + and 2 + channels and their ratio 
for j3 values in the broken phase near the weakly first order transition at L=48. The 
statistics of each simulation is also given. 



p 


m + 


m 2 + 


m 2 +/m + 


statistics 


0.5508 


0.1556(36) 


0.381(17) 


2.45(17) 


300k 


0.550875 


0.1565(56) 


0.384(16) 


2.45(19) 


200k 


0.551 


0.1837(59) 


0.444(18) 


2.42(18) 


100k 


0.552 


0.2375(42) 


0.533(36) 


2.24(19) 


100k 


0.553 


0.2900(34) 


0.660(27) 


2.28(12) 


200k 


0.554 


0.3258(60) 


0.691(57) 


2.12(21) 


100k 


0.555 


0.3502(67) 


0.847(29) 


2.42(13) 


200k 


0.556 


0.3996(85) 


0.891(35) 


2.23(14) 


200k 


0.56 


0.4965(83) 


1.204(30) 


2.43(10) 


200k 


0.562 


0.537(11) 






200k 


0.565 


0.6324(85) 






200k 


0.57 


0.702(12) 






200k 


0.575 


0.8381(78) 






200k 


0.58 


0.9358(97) 






200k 


0.60 


1.170(20) 






200k 



Since for first order transitions the critical exponent of the correlation length is v = 
1/3 [25], the scaling law 

/ Pi-Pt Y ™o+(Pi) (17) 
\P2-Pt) ~ rn 0+ ((3 2 ) { } 

must hold in the scaling region. Here m +(/3i) and mo+(P 2 ) are the fundamental masses in 
the + channel at Pi and p 2 , respectively. Comparing the scaling function with the ratios 
determined from simulations, we can estimate how large is the scaling region. Fig. [TT1 
shows that the scaling region extends at least up to (3 = 0.562 ({T t — T)/T t ~ 0.02). 

We have determined the ratio m 2 +/m + for several (3 values in the region [/3 t , 0.56] ; 
the results are presented in Table [1] and plotted in Fig. [12j This ratio remains practically 
constant in the considered region, thus suggesting that the correlation length £ 2 associated 
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to the channel 2 + (£ 2 = I/7712+) scales in the same way of the fundamental one (£q = 
l/m +). We can take as our estimation of the mass ratio the value 

— = 2.43(10) , (18) 

m + 

determined, as discussed above, by taking value and error of the point with the smallest 
error belonging to the plateau {(3 = 0.56, see Fig. H2]) . A fit of the data with a constant 
gives for this ratio the value 2.353(49) with a x 2 /d.o.f.= 0.54. 

The result given in (IT8"1) is to be taken as a prediction for the corresponding ratio of 
the screening masses in the broken phase of the (3+l)d SU(3) pure gauge theory at finite 
temperature. 



5 Conclusions and outlook 

In this work we have studied massive excitations of the 3d 3-state Potts model near 
the Ising critical point on the temperature - magnetic field phase diagram and near the 
transition point at zero magnetic field. 

We have found evidence that the mass ratio m 2+ /m + near the Ising critical point 
is compatible with the prediction from universality, thus supporting the conjecture of 
universal spectrum. 

In the broken phase of the scaling region near the transition in absence of the external 
source, we have found m2+/ m o+=2. 43(10). This result is to be compared with the cor- 
responding ratio between screening masses of (3+l)d SU(3) pure gauge theory at finite 
temperature in the broken phase near the deconfinement temperature [22] . An agreement 
between these determinations would represent an indication that universality of the mass 
spectrum holds also for weakly first order transitions. 

The present analysis could be extended with the numerical determination of higher 
masses in the + and 2 + channels and with masses in other channels, in order to carry 
out a more systematic study of universality effects. Such an extension would be certainly 
well-motivated in case of a positive answer of the comparison with the SU(3) gauge theory. 

References 

[1] A.M. Polyakov, Phys. Lett. B 72 (1978) 477. 

[2] L. Susskind, Phys. Rev. D 20 (1979) 2610. 

[3] B. Svetitsky and L.G. Yaffe, Nucl. Phys. B 210 (1982) 423. 

[4] F. Gliozzi and P. Provero, Phys. Rev. D 56 (1997) 1131. 

[5] R. Fiore, F. Gliozzi and P. Provero, Phys. Rev. D 58 (1998) 114502. 

[6] J. Engels and T. Scheideler, Nucl. Phys. B 539 (1999) 557. 



10 



[7] S. Fortunato, F. Karsch, P. Petreczky and H. Satz, Nucl. Phys. (Proc. Suppl.) 94 
(2001) 398. 

[8] R. Fiore, A. Papa and P. Provero, Nucl. Phys. (Proc. Suppl.) 106 (2002) 486. 

[9] R. Fiore, A. Papa and P. Provero, Phys. Rev. D 63 (2001) 117503. 

[10] A. Papa and C. Vena, Int. J. Mod. Phys. A 19 (2004) 3209. 

[11] A. Pelissetto and E. Vicari, Phys. Rept. 368 (2002) 549. 

[12] M. Caselle, M. Hasenbusch and P. Provero, Nucl. Phys. B 556 (1999) 575. 

[13] M. Caselle, M. Hasenbusch, P. Provero and K. Zarembo, Nucl. Phys. B 623 (2002) 
474. 

[14] R. Fiore, A. Papa and P. Provero, Nucl. Phys. (Proc. Suppl.) 119 (2003) 490. 

[15] R. Fiore, A. Papa and P. Provero, Phys. Rev. D 67 (2003) 114508. 

[16] F. Karsch and S. Stickan, Phys. Lett. B 488 (2000) 319. 

[17] H.W.J. Blote and R.H. Swendsen, Phys. Rev. Lett. 43 (1979) 779. 

[18] W. Janke and R. Villanova, Nucl. Phys. B 489 (1997) 679 and references therein. 

[19] R.V. Gavai, F. Karsch and B. Petersson, Nucl. Phys. B 322 (1989) 738. 

[20] A. S. Kronfeld, Nucl. Phys. (Proc. Suppl.) 17 (1990) 313. 

[21] M. Liischer and U. Wolff, Nucl. Phys. B 339 (1990) 222. 

[22] R. Falcone, R. Fiore, M. Gravina and A. Papa, work in progress. 

[23] R. H. Swendsen and J. S. Wang, Phys. Rev. Lett. 58 (1987) 86. 

[24] P. W. Kasteleyn and C. M. Fortuin, J. Phys. Soc. Jpn. Suppl. 26 (1969) 11. 

[25] M. E. Fisher and A. N. Berker, Phys. Rev. B 26 (1982) 2507. 



11 




30000 
25000 
20000 
15000 
10000 
5000 




1 1 1 1 1 1 


x=0.37238 - 








1,1. 


\! 


i.J 



1.5 



2.5 

M 



3.5 



30000 




Figure 3: Distributions of the order parameter M near P c . The simulations have been 
done at the points (£ c , r) with r = 0.37182 = r c , r=0.37238 and r=0.37248. These points 
lie on the first order line starting from P c . The statistics is 500k in all cases. 
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Figure 4: Effective mass in the + channel as a function of the separation between walls 
on the (y, z) plane at £ = £ c and r=0. 37248, determined from the configurations belonging 
to the "right-peak" in the thermal equilibrium ensemble. 
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Figure 5: Effective mass in the 2 + channel as a function of the separation between walls 
on the (y, z) plane at £ = £ c and r=0. 37248, determined from the configurations belonging 
to the "right-peak" in the thermal equilibrium ensemble. 
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Figure 6: Scatter plot of the complex variable S at the transition point /? f =0. 550538 [19 
at zero magnetic field. Both the symmetric and the broken phases are present. 
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Figure 7: (Top) Typical scatter plot of the complex variable S for (3 larger than 0.5508 
and zero magnetic field on a lattice 48 3 . There are no states in the symmetric phase, but 
tunneling survives between the three broken minima. 

(Bottom) Same as (Top) with the tunneling between broken minima removed by the 
"rotation" to the real phase. 
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Figure 8: Effective mass in the + channel as a function of the separation between walls 
on the (y, z) plane at (5 — 0.554 and h — 0. 
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Figure 9: Effective mass in the 2 + channel as a function of the separation between walls 
on the (y, z) plane at (5 — 0.554 and h — 0. 
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Figure 10: Fundamental masses in the + and 2 + channels as functions of /3, in the broken 
phase near T t (vertical line). 
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Figure 11: Comparison between the scaling function [(/3i — /3 t )/(/3 — j3 t )} 1 ^ and the mass 
ratio m + (Pi) / m + (P) for varying f3, with / 9 1 =0. 550875. 



19 




Figure 12: m 2 + (j3) / m + (j3) for f3 varying in the scaling region. 
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